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- - -, Abstract 

(«— ^ ' The hadron-energy evolution (Cohins and Soper) equation for ah the leading-twist transverse- 

^^ ■ momentum and spin dependent parton distributions is derived in the impact parameter space. 
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Based on this equation, we present a resummation formulas for the spin dependent structure 
functions of the semi-inclusive deep inelastic scattering. 



I. INTRODUCTION 

Semi-inclusive processes at low transverse momentum have attracted considerable inter- 
est in recent years. These processes can provide much information on aspects of perturbative 
and non-perturbative quantum chromodynamics (QCD), and on the internal structure of 
the nucleon in particular. Theoretical study of these processes began with the classical 
work of Collins and Soper in which a nearly back-to-back hadron pair is produced in e^e~ 
collision ^ij. Applications to the Drell-Yan process and the semi- inclusive deep- inelastic 
scattering (SIDIS) were made lately in [2, |3|, y]. A factorization theorem for the process 
was established Pl, involving a new class of non-perturbative hadronic observables depend- 
ing on the transverse-momentum of hadrons and/or partons: the transverse-momentum 
dependent (TMD) fragmentation functions and parton distributions. Based on the recent 
development of the gauge invariant definition of the TMD parton distributions p, |6[ , the 
proof of the factorization has been extended to the SIDIS and Drell-Yan processes at low 
transverse momentum (one the order of I^qcd) |ZI;|s|; where the TMD parton distributions 
and fragmentation functions as important ingredients were emphasized. 

On the other hand, at low transverse momentum, Pi <C Q, where Q represents some hard 
scale (e.g. the invariant mass of the virtual photon in SIDIS), the standard perturbative 
QCD (pQCD) calculations generate large logarithms a" In^" Q^Z-^l in perturbation series. 
These large logarithms must be resummed to make predictions reliable j^lSljllfll- According 
to the factorization theorem for the semi-inclusive processes, these large logarithms can be 
attributed to the hadron-energy dependence of the TMDparton distributions and fragmen- 
tation functions, which contains double logarithms [l|, |7[. The CoUins-Soper equation is 
precisely the equation which governs the energy dependence of the TMD parton distribu- 
tions and the fragmentation functions pj . The factorization proof and an extensive study for 
the energy evolution equation for the unpolarized quark distribution have been performed 
in |l|. The large logarithms mentioned above can be resummed by solving this evolution 
equation [l|, |^ . In this paper, we will follow these studies to analyze the evolution equations 
for the spin-dependent quark distributions. We will show that the original Collins-Soper 
evolution equation also applies to the spin-dependent distributions. With these evolution 
equations, we can perform the large logarithmic resummations for the spin-dependent struc- 
ture functions and asymmetries in the semi- inclusive DIS [7|, |8[ . Another important point of 
our calculations is that we use Feynman gauge, and the TMD parton distributions in this 
gauge are defined in such a way to guarantee the gauge invariance [5, 6]. In [1], a specified 



gauge (axial gauge) was used, while in [ll| the energy evolution for the Sudakov form factor 
was calculated in Feyman gauge. 

The rest of this paper is organized as follows. In Sec. II, we briefly review the factorization 
of the energy evolution equation for the TMD quark distribution following [l|, |ll| , and point 
out that this factorization works also for the spin-dependent distributions. The soft and hard 
parts in the factorization are defined. In Sec. Ill, the evolution equation for the unpolarized 
quark distribution at one-loop order is rederived. In Sec. IV, we show the results for all the 
leading-twist TMD quark distributions. In Sec.V, with these evolution equations, we show 
how to resum the large logarithms. We conclude in Sec. VI. 



II. FACTORIZATION OF THE ENERGY DERIVATIVE FOR THE TMD QUARK 
DISTRIBUTIONS 

In the non-singular gauge (e.g. Feynman gauge), the TMD quark distributions can be 
defined through the following matrix: p, lllLll2l|. 



^ rr , 0, h)Clioo; r, 0, 6±)/:.(oo; 0)V^,(0) 
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where ipg is the quark field with the Dirac- and color indices implicit. The parent hadron 
has a momentum P'^ along the ^-direction, and is polarized with a spin vector S'^. In the 
following, we use the light-cone coordinates k^ = (fc° ± k^)/\/2, and write any four- vector 
k'^ in the form of {k~, k) = {k^ , k^, k±), where k± represents two perpendicular components 
{k^,ky). The light-like vector p is chosen to be p^ = A(0, 1,0_l), and v^ is a time-like 
dimensionless (f^ > 0) four- vector with zero transverse components {v~,v^,^). We choose 
v~ ^ v^ so that v^ is very close to the light-like vector n^ = (1, 0, 0_l)/2A. The variable C^ 
is essentially the energy of the hadron, C,"^ = {2P -vY/v"^ = 2{P^Yv~ /v^ . C^ is a gauge link 
along v^, 



oo 



£^(oo; = exp I -ig I dXv ■ A{\v + I • (2) 

Here a non-light-like gauge link is introduced to regulate the li^ht-cone singularity associated 
with the plus component of the gluon momentum /"*" — >^ jjj, lll|. We avoid the use of the 
singular gauge (e.g. the light-cone gauge), because in such a gauge it is well known that 
the gauge potential does not vanish at infinity, and therefore a gauge link at infinity will be 
needed to ensure gauge invariance p. 

In the above definition, we have divided by a soft factor which is defined as 

5'(&±,/i^p) = — (0|/:t.^(6j_,-oo)/:i,j.(oo;6^)/:^jfc(oo;0)/:{)H(0;-oo)|0) , (3) 

where v is another off-light-cone vector with v^ ^ v^ and v± = 0. The parameter p is 
defined as p^ = v~v^ /v^v~ . This soft factor will also appear in the factorization theorem 
for the semi-inclusive processes at low Pj_ |7]. 

The above definition of the matrix M. is for deep inelastic scattering. For the Drell-Yan 
process, we need a different form where the gauge links are along the backward direction to 
— cxD J5|,|6|, and the soft factor has to be modified as well [8|. In the following discussions, we 
will restrict ourselves to those for SIDIS, although the same can be done for the distributions 
in the Drell-Yan process. 

The leading-twist expansion of the matrix Jvi contains eight distributions 13lll4l|. 



M 



q{x, k±) i) + TT'^9(a;, k±)a^"'k^p^ + Agi(x, k±)X'j5 ^ (4) 
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+—AqT{x, k±)-f5 ^(A;_L ■ S±) + —SgL^x, k^^Xia^^-f^p^k^i + 5qT{x, k±)icrf,^-f5p>'S] 

+ ^SqT'{x, k^)ia,,^,p^ (k^ ■ S^kl - \klS'i\ + ^qAx, k^)e^''''''^,p,k^Sp 



where M is the nucleoli mass. We have omitted the arguments C, /^ and p in the parton 
distributions on the right-hand side. The polarization vector S^ has been decomposed into 
a longitudinal component S^ and a transverse one S"^, and A is the helicity. The notations 
for the distributions follow Ref. 15| . 



The energy evolution of the TMD parton distribution is governed by the Collins-Soper 
equation [l|, 

d 

C^/(a;, &, /i, xC, p) = {K{h, p, p) + G(xC, p, p)) f{x, b, p, xC, p) , (5) 

where K and G depend on soft and hard physics, respectively. In the following we will 
show that this equation is in fact valid for all the leading-twist quark distributions in the 
impact parameter space. It is important to note that the above factorization of soft and 
hard parts for the energy derivative is true only for the leading power (1/C) contribution, 
and any higher power corrections have been neglected [l|, lll|. So, the evolution equations 
like Eq. (0) and here after are only valid in the leading power of l/C Another point we want 
to point out is the p dependence in K and G, which was absent in [l|, lll|. The p dependence 



comes from our definition of K in Eq. (|13|) below to avoid possible light-cone singularity in 
higher loop calculations. However, the p dependence of K and G cancels out, and energy 
derivative of the TMD parton distribution does not depend on p, (see detailed discussion 
below) . 

Since there is no energy dependence in the soft factor S", the only source for energy 
dependence comes from the numerator in the matrix M. in Eq. (U). We call the numerator 
the un-subtracted TMD parton distribution Q(a;, A;_l,xC) as in ^^J. So, the Collins-Soper 
equation for any un-subtracted distribution reads 

d 

C—J^{x, b, /i, xC) = {K{b, p, p) + G{xC, p, p)) Tix, b, p, xQ . (6) 

where K and G depend on soft and hard physics, respectively, as we shall explain. Since any 
un-subtracted distribution does not depend on the soft part, then there is no p dependence 
in J^. This means that sum {K{b, p, p) + G{xC, p, p)) is p-independent while each term could 
be p-dependent. 

The different TMD quark distributions can be obtained from Eq. (4) by making appro- 
priate spin projections. For example, the unpolarized quark distribution is related to, 

Q(x, fc^, /i, <) = \j ^^9^e-«-^"+^^-^^- (P |^,(r, b^)Ch^C,U^) I P) . (7) 

Since C,"^ = {2v ■ P)^/f^, using the chain rule, the derivative on ( can be related to the 
derivative on v, 

d d 

where 5v is another dimensionless vector: 5v~ = v~ , 6v^ = —v~^, and 6v± = 0. So that, 
we have Sv'^ = —v"^ < and Sv ■ v = 0. From Eq. ((Tj), we see that the only dependence 
of V comes from the gauge link Cy. So, we have the following differential equation for the 



v-k+ie^ ^ ' v-k 



v-kSv'^—Sv-kv^ 
^ (v-k+ie)'^ 



5_^M 



FIG. 1: Feynman rules for the eikonal vertex in the TMD parton distribution (a) and its derivative 
(b). 



un-subtracted TMD parton distribution, 
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The relevant Feynman rules for the eikonal vertex from the gauge link contribution are 
shown in Fig. 1 for the TMD parton distributions (a), and the derivative Cw (b). 

Similar to the structure functions for the semi- inclusive DIS discussed in l7l ^1 , the deriva- 



tive C^ on the un-subtracted distribution Q receives contributions from three different 
regions of the gluon momentum: soft, hard, and coUinear contributions, respectively. 

First, let us examine the coUinear contribution to the derivative. From the above, the 
special vertex reads, 

V ■ k5v^ — 5v ■ kv^ , ^ 

[v ■ k + le)-^ 

If the gluon momentum k is in the collinear region, i.e., k oc P, we have k^ r^ Q, k^ r^ XQ, 
k± ~ vAQ, where A is a small parameter. For such momentum, the above vertex will lead 
to 

[V ■ k + ley 

and their contribution will be suppressed if we contract the above vertex with the collinear 
momentum in the jet part of the TMD parton distributions. 

Thus, the derivative C^ on the parton distribution will receive contributions only from 
the soft and hard regions of the gluon momentum. A detailed analysis of these contributions 
leads to a factorization of the derivative of the TMD distribution as illustrated in the second 
line in Fig. 2 |l|,lll|. This can be represented by the following differential equation. 



C-^Q{x,k^,xO = j[K + G]0Qix,k^,xC) 



(11) 
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FIG. 2: Factorization for the derivative d/d In C, on the un- subtracted TMD parton distribution. 

where the soft part is called K, and hard part G; and the label ® means the momentum 
space convolution. Again, we emphasize that the above factorization is only valid for the 
leading power contribution. After applying a Fourier transformation to the coordinate space, 
the above differential equation reads. 



d 

C^Q(s, 6, xQ = [K{b, /i, p) + G(xC, /i, p)] X Q{x, b, xQ 



(12) 



where the convolution in momentum space becomes products in the impact parameter b- 
space. The soft-part K depends on b, while the hard part G depends on hard scale (; and 
both of them depend on the renormalization scale p and p, but the sum does not. The 
above equation is valid for any value of b. If b is small as 1/6 ^ ^qcd, we can further have 
a factorization for the TMD parton distribution which depends on the integrated parton 
distributions, and then we can get another form for the evolution equation with an extra 
term pj. 

The soft part can be calculated from the Feynman diagrams by using the Grammer- 
Yennie approximation [l6|. In |7[, we demonstrated how to get the soft contribution in the 
TMD parton distribution. Here, we follow the same procedure, and define the soft part for 
the derivative of the TMD quark distributions by the following matrix element. 



Kib,p,p) = (0 



ig) / d\[5vA{\v) + \5vdA{\v)-v] 



/:t(6j_,-cx))£|;(oo,6j 



+ {ig) / d\ [Sv ■ A{\v + b^) + \5v ■ dA{\v + b±) ■ v] \ C^{oo, 0)C^{0, -oo) 







Clib±, -oo)Cl{oo, b^)Cy{oo, 0)/:s(0, -oo' 









(13) 



where we have introduced the off-light-cone vector v defined above to regulate possible 
light-cone singularities [7[. Notice that this definition is different from that in [ll| where a 
light-cone vector p = (0, 1, 0_l) was used instead of our off-light-cone vector v. Introducing v 
leads to the p dependence of K. Since the leading order (one-loop) calculation has no light- 
cone singularity, there is no difference between using v and p at this order, and the result 
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does not depend on p [11]. However, at higher- loop, there might be hght-cone singularity 
associated with the gluon momentum /"*"—!> 0. We note that so far there is no explicit 
calculation of the soft factor K beyond one-loop order in non-abelian gauge theory (a QED 
result has been given in [ll|) . It is not clear at present that we can take the light-cone 
limit for higher loop calculations. So, it is necessary to include the off-light- cone vector v in 
the formal definition. If there is no light-cone singularity, we can take the light-cone limit 
({)—>• p and p — > oo). Like in the factorization of SIDIS structure function [7|, p is just a 
parameter which separates the hard and soft physics, and it does not affect any prediction 
power in the resummation formalism. Normally, we should take p ^ 1, while in practice 
we can choose p between 3 and 10 if there is any p dependence in K to avoid the large 
logarithms associate with p. 

Comparing Eq. (fT^ with Eq. Q, we find that the soft factor in the derivative can be 
related to that in the parton distribution, 

^c- "• "> = w^f-'t"-"' ^' "> = ^fi''-"' '' "> ■ ("' 

Our one- loop result below verifies this relation. Since the distribution is p-independent, the 
p dependence in K must cancel the p dependence in G{x(, p, p). 

The renormalization scale dependence of the soft-factor K is determined by the cusp 
anomalous dimension |l7| . 

^Tp^ = -^^ ' ^'^^ 

which is a series in a^ and free of infrared singularities. The hard part G can be calculated 
through a systematic subtraction, as will be illustrated by the one-loop example in the next 
section. From the definition, it is obvious that the soft factor is spin-independent. 

The above analysis of the factorization of energy derivative and the definitions of the hard 
and soft parts can be extended to the spin-dependent TMD quark distributions, because they 
all come from the same matrix A4 Eq. (1) with different spin projections, and the arguments 
supporting the factorization can be generalized to all the leading-twist distributions. We will 
explicitly show this in more detail in the following one-loop calculations, and give general 
argument for all orders. 



III. RE-DERIVATION OF COLLINS-SOPER EQUATION FOR UNPOLARIZED 
QUARK DISTRIBUTION 

In this section, we demonstrate how to calculate, in the present formulation, the Collins- 
Soper evolution kernels K and G for the unpolarized quark distribution Q at one-loop order. 
In next section, we will discuss the spin dependence and present the evolution equations for 
all the leading-twist TMD quark distributions. 

We first calculate the leading contribution to the soft part K, linearly proportional to the 
strong coupling constant a^. The relevant diagrams are shown in Fig. 3. Fig. 3(a) vanishes 
because 6v ■ v = 0, and 3(c) vanishes for the same reason and for being sub- leading order in 
1/C^. The contribution from Fig. 3(b) is 

a Gp Li 
K{b,p,p)\{ig,3(^b) = ^In— , (16) 



+ 




+ 
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(a) 



(b) 



(c) 




FIG. 3: One-loop diagrams for the soft part K . The mirror diagrams are not shown, hut included 
in the result. 




FIG. 4: One-loop diagrams for the Collins-Soper evolution for the un-suhtracted TMD parton 
distributions. 

where we have included a factor of two to account for the two vertex correction diagrams. 
A is the gluon mass, introduced to regulate the infrared singularity for individual diagrams. 
However, the sum of all contributions is free of the infrared divergence. Fig. 3(d), together 
with its mirror diagram, contributes 



Summing up, we get. 



^(&,/i,p)|fig.3(d) 



K{b,fx,p) 



a.Cf 



TT 



In 



4 



62 A2 



27, 



a.C 



s^F 



71 



In ^ + 2jE 



(17) 



(18) 



which agrees with the previous calculations [l|, |ll|, |l8[ . Moreover, up to first order it has no 

dependence on p, which means that we can take the light-cone limit for v (i.e., let v = p) at 
this order. We remark here that higher order calculation of K{b, p, p) may show explicit p 
dependence. Comparing this result with the soft factor S{b, p, p) at one-loop order calculated 
in |2|, Eq. (fT^ is clearly satisfied. The one- loop cusp anomalous dimension is. 



d 



a.Cf 



-IK 



TT 



(19) 



which is well known. As we stated in the previous section, the soft part is spin independent, 
hence the above result for K is the same for all the leading order TMD quark distributions. 
We now calculate the complete one- loop contribution to the right-hand side of the Collins- 
Soper equation, from which we will subtract the above soft contribution to get the hard part 
G. All the one-loop diagrams are shown in Fig. 4. The contributions from Fig. 4(a) and (c) 
vanish because of the same reason as that for Fig. 3(a) and (c). 
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The contribution from Fig. 4(b) in momentum space reads, 

^Q(x,fc^,xC)|fig.4(b) = ^ (l-ln^^ Q{x,k^,xC) , (20) 

which is the same as that in the impact parameter b space. The exphcit dependence on 
C and A indicates the presence of both hard and soft contributions. Subtracting the soft 
contribution in Eq. (fTT)|) . we get the hard part G as, 

G(xC,/i) = ^(l-ln^) , (21) 

which depends on the hard scale ( and the renormahzation scale /j,. 

Fig. 4(d) is dominated by the contribution from the soft gluon momentum region: q~^ <^ 
k~^, where q is the gluon momentum and k is the quark momentum. After the soft-gluon 
approximation, we get 

•9 _^ , ... dsCp f (Pq± 



Q{x,k^,xC)\&gA{d) = ^^ / r^—^Q{x,k^-q^,xC) ■ (22) 



91nC 
Fourier-transforming to the b space, we have 

^^Q(x,6,xC)|fig.4(d) = —^ 



'^¥x^-^^^ 



Q{x,b,xC), (23) 



which can be reproduced by the soft factor K from Eq. (fT7|) (Fig. 3(d)). 

The above results show that the factorization is valid at one-loop order with the sum of 
K and G reads, 

K{b, /.) + GixC, /x) = -^^ In ^e^^--i . (24) 

This result agrees also with that in [7] , where the distribution itself was calculated to one- 
loop order. 



IV. COLLINS-SOPER EVOLUTION EQUATIONS FOR SPIN-DEPENDENT 
TMD DISTRIBUTIONS 

In this section, we study the energy evolution of the spin-dependent TMD distributions. 
Except for the unpolarized q{x, k±), all other leading-twist TMD quark distributions depend 
on the polarization of either the initial hadron or the probing quark. The polarized distri- 
butions can be obtained from spin projections of the matrix Ai Eq. (^. There are three 
leading-twist projections. 



7 
7^7'7' 



qix,k^), qTix,k±) ; 

AqL{x,k±), AqT{x,k±) ; 

SqT{x,k_i_), 6qL{x,k_i_), 6q{x,k_L), 6q'rp{x,k±) , (25) 



corresponding to the unpolarized, longitudinally-polarized, and transversely-polarized quark 
distributions, respectively. Moreover, different distributions may have different k_j_ depen- 
dence. For example, so-called k±-even (under the exchange k± — > —k±) quark distributions, 



q{x,kj_), AqL{x,k±), and Sqt(x , k\ ), correspond to the unpolarized, helicity, and transver- 
sity distributions, respectively [19J. These distributions survive after integrating over k±. 
The other five distributions are associated with the k±-odd structures, and vanish when k± 
are integrated over. 

In the following, we use the Sivers function ^C| as an example to demonstrate the calcu- 
lation of the energy evolution kernel. The results for the other distributions can be obtained 
similarly. The Sivers function results through 7+ projection from the matrix A4 Eq. (^, 



QT{x,k±,fj,,x() 



2e'^Sikj 
X (PSi 



(27r)3 

^jrA)ch''c,M^) PS, 



spin dependent part 



(26) 



where the explicit transverse momentum and spin dependence has been included. The Feyn- 
man diagrams for the energy derivative are the same as those for the unpolarized distribution 
discussed in last section; and Figs. 4(a) and (c) vanish as before. The contribution from 
Fig. 4(d) reads. 



d 



ainC 



e'^SikjQ^{x, fc±,a;C)|fig.4(d) 



a.C 



s'-^F 



27r2 



(fa 
^ \ e'^Si{k - q)jQj,{x, k^ - gl, xQ 



(27) 



where again we have made the soft approximation. To find the above result, we have applied 
the 7"*" projection to the quark matrix Ai. Moreover, because the Sivers function is spin- 
dependent and associated with a A;_L-odd structure, only such structure is isolated and kept. 
Fourier transforming to the impact parameter space, we get. 



a6Qr(a;,o,a;C)|fig.4(d) = 



ainC 



TT 



In 



62A2 



-27b 



dlQT{x,h,xQ 



(28) 



where dl = d/dbi is a derivative on the Sivers function. Apart from the explicit derivative, 
the above contribution is the same as that for the unpolarized distribution in the previous 
section. 

The contribution from Fig. 4(b) is also the same as that for the unpolarized quark dis- 
tribution, 



d 



ainC 



QT{x,k±, xQlfigAih) 



""'^^ ^l-\n^]QT{x,k^,xC) 



TT 



A2 



(29) 



and similar equation holds in the 6-space. Combining the above results, we have the entire 
energy evolution of the Sivers function at one-loop order. 



d ■ 
CgjdlqTix, b, /i, xC, p) = {K{b, /i, p) + G{xC, /i, p)) dlqrix, b, /i, xC, p) 



(30) 



where K and G are the same as those for the unpolarized distribution. 

The above analysis can be repeated for all other leading-twist quark distributions, and the 
one-loop evolution kernels are found again to be the same as those for the unpolarized quark 
distribution. There are two important features supporting the above finding. First, the 
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leading-twist projection matrices 7+, 7^7^, and 7+7*7^ lead to an identical trace. Second, 
there is no mixing between the different leading-twist distributions. This latter property, 
however, does not hold for the evolution of the higher-twist distributions, a topic beyond 
the scope of this paper. 

In fact, we argue that, to all orders in perturbation theory, the same evolution kernels 
determine the energy evolution of all the leading-twist distributions. First of all, from the 
factorization of the energy derivative discussed in the Section II, only soft and hard regions 
contribute to the evolution. For the soft part, there is no spin dependence, as is clear 
from its definition Eq. (fT^ . Therefore, for any leading-twist distribution, the soft part 
of the evolution kernel is the same. Second, the contribution from the hard part is also 
spin-independent, because the hard contribution is calculable in perturbative QCD and the 
perturbative processes for massless quarks conserve helicity. Any spin projection will lead 
to the same Dirac algebra if there is no mixing between different distributions. 

To summarize, the CoUins-Soper evolution kernel has no spin dependence for the leading- 
twist TMD quark distributions. For /c^-even istributions, the following evolution equation, 

d 
C^/(a;, b, /i, xC, p) = {K{h, /i, p) + G(xC, /i, p)) f{x, b, p, xC, p) , (31) 

holds for f = q, Ag^, and 6qT- For k±-odd quark distributions, we have 

d ■ 
C-Q^dlf{x, 6, /i, xC, p) = iK{b, /i, p) + G{xC, p, p)) dlf{x, b, p, xC, p) , (32) 

which works for / = qx, Ag-p, 6q, and 6qL. Finally, for 6q'rp, 
d 



C^ [dldi - 6'^dl/2j 6q'Ax, 6, p, xC, p) 

= {K{b, p, p) + G(xC, p. p)) (dldi - 5''dll2) Sq'Ax, b, p, xC, p) . (33) 

As a final remark, we would like to point out that the energy evolution equations for the 
TMD quark distributions in the Drell-Yan process will be the same as the above for the DIS 
process. This is because we have the universality for the parton distributions between the 
two processes p, |6|. 



V. RESUMMATION FOR THE STRUCTURE FUNCTIONS IN POLARIZED 
SEMI-INCLUSIVE DIS 

In the physical cross sections with two widely separated scales, say. Pi and Q, there are 
large double logarithms of the type a" In "^ Q'^IP'\_ as well as sub-leading ones. To have a 
reliable theoretical prediction, one has to re-sum these contributions. In this section, we 
perform the resummation for the large logarithms in polarized semi-inclusive DIS by solving 
the CoUins-Soper evolution equations obtained in the last section. 

In J8| , we have obtained the factorization formulas for the various structure functions. In 
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the impact parameter b space, they can be expressed as 



Flj'lib) = q{xB.Zr,h)q{zuMS^{h)H^u}jm , 

pLLih) = AqL{xB,Zhb)q{zH,b)S+ib)HLL{Q' 



^(1) 



F^m 



^(^)i 



Mzh 
-1 



-dl 



{dlqrixB, Zhb)) q{zH, b) S+ {b) H\;},{Q^) 



F^^m ^ j^^dl 



6qTixB,Zhb) {dl5qizj,,b)) S+{b)H\j%{Q 



Fijiib) 
FuL{b) 
Fijkb) 



Mzh 



dl [{dl/\qT{xB,Zhb)) q{zH,b)S+{b)HLT{Q^)] 



dldi - dl5'^ 



-'b'^b 



MMhZh 



{di6q{xB,Zhb)) {diSq{z,,b)) S+{b)H\flj{Q 



" ' " [{diSqLixB, z,b)) {di6qiz,, b)) S-^ib)HuLiQ' 



MMhZh 



Mm^zi 

{dt5q{znM)S^{b)HSl{Q^ 



{dldi-6^^d^j2)6q'j^{xB,z,b) 



(34) 



where we followed the notations used in J8| , and the parton distributions and fragmentation 
functions are calculated at the energy scale x^C^ = C^/^I = pQ^- F{b) are the Fourier 
transformation of the structure functions in the impact parameter space. For the unpolarized 
structure function, we define 



^(1), 



(1), 



F^-'ib) = d'P,^F^'>{P,^)e 



JPh I b I 



(35) 



and similarly for Fll. Others denote transverse- momentum- weighted Fourier transforma- 
tions, for example, 

F«(6) = [ d'PH^\Ph±\Fij'kPk±)e''''^^-'^ , (36) 



^(2) 



(2) 



and similarly for F^-f, and Fit- For F^^j and Ful, we define 



For Fjjj^, we define 



Flj^ib) 



Fl?kb) 



d'P.AP ^'^^'^ 



h±\-rh±\ 



F^Z{Pk±)e 



iPh I b I 



- I d'P,^\P,^\'FljkPh±)e''''^^'' 



(37) 



(38) 



In [7|, the large logarithms have been re-summed for the unpolarized structure function 
by solving the relevant Collins-Soper equation for the TMD quark distribution and frag- 
mentation function. From Eq. (|34j) and the results of the previous section concerning the 
Collins-Soper evolution for the polarized quark distributions, we conclude that the polarized 
structure functions have the same evolution equation as the unpolarized case. In the follow- 
ing, we take the structure function -F^j., depending on the transversely-polarized nucleon 
spin, as an example to demonstrate the re-summation procedure. 
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Rewriting the structure function FJ)^ as, 

Fi^kb, Q') = -W^J^l^UTih. Q') ' (39) 

where 

:FuT{h, Q^) = {dlgrixB, Zhb)) q{zh. h)S+{h)H^^l{Q^) . (40) 

Since the Sivers function qt and the unpolarized fragmentation function q obey the same 
Colhns-Soper evolution equation as that for the unpolarized quark distribution, we have the 
following derivative equation respect to Q^ |2,l2|, 

d 
Q^-k7^^ut{xb, Zh, b, Q'^) = [Kihji, g{fi),p) + G'^j,{Q/fi, g{fi),p)] J^ut{xb, Zh, b, Q^) , (41) 

where K is the same as before, and G'jjrp contains additional contribution from the hard 
part. The sum of K and G'^rp has no dependence on p although they separately might have. 
The solution to this differential equation has the following form J2|, 

J^uAxB, ZH, b, Q2) = J^ut{xb, Zh, b, pl/CDe-^^"^"'''^'^ , (42) 

where the distribution and fragmentation function are evaluated at x^C^ = C'^/z^ = p/xi/C|. 
The Sudakov suppression form factor reads. 



S{Q\plb,C2)= f 

J a 






\n[^]A{bpL,P,p) + B{C2,bpL,p,p) . (43) 



Here C2 is a parameter in the order of 1, and pi is a low-energy, but still perturbative, 
scale. The functions A and B have perturbative expansions in a^, A = X]n^''"H'^s/^)" ^^^ 
B = Y.nB^''\(^s/T^T- They are defined as 

d 
A{bpL,p,p) = jk{p) + P-g-K{bpL,g{p),p) , 

BiC2,bpL,p,p) = -2K{bpL,g{p),p)~2G'jjT{l/C2,g{p),p) . (44) 

The y4-function is the same as that for the unpolarized case J2, 111 since it comes only from the 
spin-independent soft part. On the other hand, the 5-function contains contribution from 
the unknown hard part HJjj^, and hence it could be different from that of the unpolarized 
one. 

Substituting the above results into Eq. (|H^. we will get |7[ 

Fij'kb, Q') = -^f, \{cf,qT{xB, Zhb)) qiz,, b)S^{b)Hi;kA/C!)e~^^^'''^l'''^] . (45) 

From the above derivations, we confirmed that the p dependence in the soft factor K does not 
affect the resummation of the large logarithms. This is because the resummation concerns 
the logarithms of the form In^ Q^b^ in impact parameter space, while p is just a parameter 
separating the hard and soft physics. In addition, there is no p dependence in the Sudakov 
suppression form factor S. 
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The re-sumination for all the other structure functions can be done in a similar way. With 
these re-summation formulas, one can further study the Q^ dependence of the P± spectrum 
of the polarized cross section and asymmetries in the semi- inclusive DIS J2l|. 

If Q^ is not too large, for example, at the order of tens GeV^, it is legitimate to re-sum 
only the leading double logarithms (DL) [2, Isl, llOf • In this approximation, we only need to 
take into account the first term in the expansion of A function, neglecting the contribution 
from B. Since A^^^ = 4/3 from the cusp anomalous dimension 'Jk, the Sudakov suppression 
factor reduces to 

5(-)(g^/.i,C,) = l^^lnf^) , (46) 

in the DL approximation. It only depends on Q^ and /i|, but not on b. Moreover, since the 
the cusp anomalous dimension is the same for all the leading-twist TMD quark distribution, 
the Sudakov suppression factor will be the same for all the leading-twist polarized structure 
functions in (J2D). One of the consequences is that one can predict the Ph± spectrum at 
higher Q^ for semi-inclusive hadron production in DIS from the result at lower Q^. For 
example, in the DL approximation, the unpolarized structure function Fut and the single 
transversely polarized structure function FJj^ have the following dependence on Q^, 

Fuu{xB,Zh,PH±,Q') = Fuu{xB,z,,P,^,fil)e-^"'''^'^''^^'^^ 

The polarization asymmetry (the ratio of these two structure functions) as a function of 
Ph± will remain the same for different Q^ at fixed xb and z^. This constancy in Q^ has 
been seen from the comparison of the HERMES data at HERA with that of CLAS at JLab 
on various spin asymmetries, with the average Q^ varying by a factor of three J22]. The 
above analysis applies to all the leading-twist polarized structure functions and polarization 
asymmetries. It will be interested to test this prediction based on DL approximation with 
future DIS experiments at different Q^. 

We notice that the above formalism also applies for the Drell-Yan process, which have 
plenty data at low transverse momentum and not very high Q^ J23|. It will be useful to 
compare the above DL approximation prediction with these experimental data, and gain 
insight for the transverse momentum dependence for the TMD quark distributions. We 
will carry this out in a futu re p ublication. This approach is different from what has been 
done so far in the literature 2J, |25|, |26|, |27|, |23| , where the predictions solely depend on the 



integrated parton distributions and the CoUins-Soper-Sterman formalism j^, and one needs 
nonperturbative parametrization for the large b behavior of the distributions and the evolu- 
tion as well. In our approach, the TMD parton distributions are important nonperturbative 
ingredients. 

If Q^ is very large (e.g., for W and Z bosons production ), the above approximation 
breaks down. One has to take into account sub-leading logarithmic contributions, perhaps 
up to v4(2) and B^^^ in the expansion of the functions A and 5 [i H H H H H Q . 



VI. CONCLUSIONS 

In this paper, we have studied the CoUins-Soper energy evolution equation for all the 
leading-twist TMD quark distributions. The evolution equation has contributions from both 
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the hard and soft regions of the gluon momentum. Since both parts are independent of the 
quark hehcity, the evolution kernel is spin independent. Based on the evolution equation, 
we can perform re-summation for the large double logarithms in the polarized structure 
functions. 

We thank H. Avakian for the communications about the transverse momentum spectrum 
of the polarized asymmetries. A. I., X. J. and F. Y. were supported by the U. S. Department 
of Energy via grant DE-FG02-93ER-40762. J.P.M. was supported by National Natural 
Science Foundation of P.R. China through grand No. 19925520. X. J. is also supported by 
an Overseas Outstanding Young Chinese Scientist grant from NSFC. 
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